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I.  INTRODUCTION 


The  equations  describing  radar  detection  probabilities  of  constant 
amplitude,  Rayleigh,  and  dominant  plus  Rayleigh  targets  have  been  well 
known  for  several  decades.  For  these  scatterers,  plots  showing 
detection  probability  versus  signal-to-noise  ratio  for  different  values 
of  false  alarm  probability  and  number  of  pulses  integrated  have  been 
documented  by  many  investigators. 

Woodward  [1],  Swerling  (2),  and  Marcum  [3],  [4]  introduced  these 
basic  stochastic  models  for  radar  scatterers.  These  early  works  are 
presented  along  with  plotted  data  in  an  excellent  treatment  of  the 
subject  by  DiFranco  and  Rubin  [5] .  Many  ramifications  pertaining  to 
system  modes  and  performance  have  been  investigated.  Coherent  and 
incoherent  integration  and  ramifications  thereof  have  been  treated  in 
many  papers  [6]  -  [10].  These  studies  have  contributed  much 
understanding  on  such  issues  as  fluctuation  loss,  collapsing  loss,  and 
diversity  gain. 

North  [6]  introduced  the  concepts  of  matched  filters,  effects  of 
post-detection  integration,  statistical  criteria  for  detection, 
optimality  of  a  square-law  detector,  and  fundamental  significance  of 
signal-to-noise  energy  ratio.  Blake  in  [7]  presents  comparisons  between 
video  integration  of  several  pulses  and  cumulative  detection  based  on 
independent  threshold  crossings  on  successive  scans.  Miller  and 
Bernstein  [8]  treated  signal-to-noise  ratio  requirements  for  coherent 
integration  along  with  detection  probability  where  a  number  of  parallel 


filters  were  used  to  cover  a  given  band  with  the  overall  false  alarm 
probability  held  constant.  Many  other  works  on  the  effects  of  frequency 
diversity  may  be  found  in  [9].  Recent  work  by  Dana  and  Moraitis  [10] 
included  the  probability  of  detecting  a  Swerling  I  target  on  two 
correlated  observations.  A  second  dwell  was  used  to  verify  the  initial 
detection  and  minimize  false  alarms.  This  technique  can  be  applied  to 
the  radar  that  operates  with  an  electronically  steered  beam, 
range-ambiguous  waveforms,  and  frequency  agility.  The  conditional 
probability  of  a  second  dwell  detection  given  a  first  dwell  detection 
(and  the  probability  of  at  least  one  dwell  detection  out  of  two  dwells) 
was  compared  with  the  probability  of  independent  dwells  when  two 
separate  transmission  frequencies  were  used. 

Many  sequential  detection  approaches  involving  integration,  either 
coherent  or  incoherent  integration,  have  been  investigated  for  multiple 
range  and  beam  positions  [11]  -  [13].  These  approaches  assumed  that,  at 
most,  one  target  is  present  with  some  weighted  probability  over  the 
sector  being  investigated.  For  example,  Finn,  [11]  conducted  a  study 
dealing  with  point  targets  on  two-step  and  three-step  sequential 
detection  systems  using  different  thresholds  and  energy  levels. 

Preston  [14]  made  comparisons  with  nonsequential  performance  to  show 
improvement  as  a  function  of  the  number  of  range  cells  when  the 
sequential  process  is  optimized  for  multiple  range  cells. 

Many  studies  have  treated  the  over  resolved  effects  of  targets. 
Nitzberg  [15]  evaluated  detection  performance  of  a  single  statistic 


output  derived  from  a  broadband  waveform  incoherently  integrated  over 
each  range  cell  spanning  the  target.  Also  evaluated  was  the  detection 
performance  for  incoherently  integrated  frequency  diversity  waveforms 
having  the  same  bandwidth  occupancy  as  the  broadband  waveform.  It  was 
shown  that  at  the  minimum  signal-to-noise  ratio, ^  detection 
performance  differed  by  only  0.8  dB  for  the  two  approaches.  Recently, 
Hughes  [16]  presented  a  high-resolution  radar  detection  strategy  in 
conjunction  with  a  constant  false  alarm  rate  (CFAR)  processor.  Hughes 
details  a  comparison  of  single  range  cell  detection  with  incoherent 
integrated  contiguous  range  cell  detection.  Integrated  range  cell 
detection  provided  improved  detection  capability  except  when  the  target 
is  composed  of  a  single  strong  specular  point. 

In  this  paper,  a  detection  is  based  on  at  least  one  threshold 
crossing  out  of  the  M  range  cells  spanning  the  target;  likewise  for  a 
false  alarm  when  the  target  is  absent.  Sequential  detection  or 
incoherent  range  cell  integration  was  not  considered  here.  Neither  were 
there  any  procedures  for  correlating  the  over  resolved  target  signature 
for  maximizing  system1 output.  Probability  of  detection  as  a  function  of 
signal-to-noise  ratio  is  provided,  with  the  number  of  over  resolved 
range  cells,  M,  as  a  parameter.  False  alarm  probability  for  the  spanned 
target  is  held  constant.  We  shall  proceed  by  describing  our  approach. 


1 

The  signal-to-noise  ratio  required  for  a  certain  level  of  detection 
performance  when  the  waveform  bandwidth  just  resolves  the  individual 
subscatterers. 
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II.  APPROACH 


First  the  threshold  value  is  determined  from  a  selected  false  alarm 
probability.  This  involves  the  use  of  analytical,  numerical,  or 
graphical  false  alarm  and  threshold  relationships.  The  various 
detection  probability  curves  as  a  function  of  signal-to-noise  ratio  can 
then  be  plotted.  Each  set  of  curves  represents  the  detection 
probability  of  a  particular  scatterer  with  a  specified  false  alarm 
probability  and  resolution  index,  M  *  1,  2,  5,  and  10,  as  a  parameter. 
When  the  target  occupies  one  range  cell  (M  =  1)  the  target  is  said  to  be 
under  or  properly  resolved.  When  the  waveform  is  such  that  the  target 
consists  of  many  subscatterers  (M  >  1)  the  target  is  over  resolved.  M 
shall  designate  the  resolution  index. 

The  data  presented  in  this  paper  shows  the  effects  of  bandwidth  (or 
resolution)  on  the  probability  of  detection  where  bandwidth  is  related 
to  the  resolution  index,  M  *  1,  2,  5,  and  10.  In  any  one  set  of  curves, 
all  other  factors  are  held  constant  over  the  different  H  values,  i.e. 
false  alarm  probability  for  the  entire  target,  total  energy 
backscattered,  and  the  number  of  pulses. 

The  range  dependent  average  cross  section  is  assumed  to  be 
uniformly  distributed  for  all  scatterers.  Although  this  may  seem  to  be 
a  special  case,  it  does  lead  to  an  upper  bound  in  performance  (the  most 
optimum  detection  probabilities)  for  increasing  M.  Scatterers 
arbitrarily  distributed  over  their  M  range  cells  yield  less  optimistic 
detection  probabilities  for  increasing  M.  Such  behavior  is  expected  due 


to  the  introduction  of  collapsing  noise  into  the  receiver.  Collapsing 

noise  is  created  when  some  of  the  received  data  samples  have  lower 

,  2 

signal-to-noise  ratio  than  that  of  other  received  data  samples. 

Ill .  FALSE  ALARM  AND  THRESHOLD  DETERMINATION 

In  establishing  radar  detection  criteria  it  is  necessary  to  first 
determine  an. allowable  false  alarm  probability,  P^.  From  this  value 
one  can  determine  a  threshold  setting  for  the  receiver  output.  A  false 
alarm  occurs  when  the  target  is  absent  and  the  noise  at  the  receiver 
output  exceeds  the  threshold.  If  the  target  is  over  resolved,  i.e. 

M  >  1,  then  a  false  alarm  occurs  when  the  noise  from  at  least  one  of  the 
M  over  resolved  range  cells  exceeds  the  threshold.  Mathematically  the 
probability  of  such  a  false  alarm  for  the  single  pulse  case  is  written 

Pfa  *  1  •*  l  1  -  exp(  -  )  ]M  (!)• 

where  denotes  the  threshold  to  noise  energy  setting  per  cell.  The 
expression  for  the  multiple  pulse  case  of  N  incoherent  integrated  pulses 
is  written 


2 

See  reference  [17]  for  a  discussion  on  collapsing  noise. 


Pfa  "  1  "  I  « 


,  N-l  )  ]“ 


(2) 


\F 


I(.,.)  is  the  incomplete  gamma  function.  In  the  above  equations  the 
threshold  setting  will  be  dependent  on  ,  the  false  alarm  probability, 
M,  the  number  of  range  cells  into  which  the  target  has  been  over 
resolved,  and  N,  the  number  of  pulses. 


Analytical,  numerical,  and  graphical  solutions  (in  that  order 
whichever  the  case  may  be)  were  used  to  determine  the  receiver 
threshold.  Threshold  values  were  determined  for  each  of  the  given 
values  of  false  alarm  probability,  resolution  index,  M,  and  number,  N, 
of  incoherent  integrated  pulses  per  cell. 


IV.  DECISION  RULE  FOR  SCATTERER  DETECTION 

A  detection  is  declared  whenever  at  least  one  threshold  crossing 
occurs  out  of  the  M  subscatterer  range  cells  spanning  the  target.  The 
example  in  figure  1  illustrates  this  point.  For  scatterer  A,  only  one 
threshold  crossing  out  of  five  subscatterers  investigated  is  observed 
and  so  scatterer  A  is  detected.  For  scatterer  B,  three  threshold 
crossings  are  observed,  thus  scatterer  B  is  also  detected.  Each 
subscatterer  response  in  figure  1  represents  a  N  pulse  integrated 
variable  where  N  >  1 .  The  receiver  for  such  a  decision  rule  is  shown  in 


figure  2. 
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Figure  1.  Example  of  the  Threshold  Crossings  for  M  *  5  Subscatterers 
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The  index,  j  (j  =  k+1 ,  k+2,  . ..,  k+M  for  any  k) ,  denotes  the  j  th 
range  cell  while  the  index,  i  (i  =  1,  2,  3,  . N) ,  denotes  the  i  th 
pulse.  A  target  is  declared  when  at  least  one  Y(j)  exceeds  the 
predetermined  threshold  Y^,  i.e. 

M 

(J  Y(j)  >  Y^  :  Hi  Target  declared 

3=1 

M 

D  Y(j)  <.  Y^  ;  Hq  Target  not  declared 

j=l 

If  the  probability  of  a  single  threshold  crossing  out  of  M 

subscatterer  cells  is  P_u,  then  from  the  binomial  formula  the 

UN 

probability  PD  of  at  least  one  crossing  out  of  M  cells  is 

PD'1-(  1-PD«)M  (3>- 

The  computation  of  P^  depends  on  the  scatterer  model,  waveform, 
bandwidth  (resolution  index),  threshold  setting,  signal-to-noise  ratio, 
and  the  number  of  pulses  being  considered.  In  the  following  sections  we 
provide  the  expression  for  P^  under  such  conditions.  Expressions  for 
the  subscatterer  probability  of  detection,  P  are  obtained  from  [5] 
with  appropriate  modifications  on  signal  energy,  E.  The  total  average 
energy  backscattered  from  all  the  subscatterers  will  be  held  constant. 
This  is  written  as  follows 


M 


(4). 


e-2  -mem 

i-1 

where  E  denotes  the  average  energy  backscattered  from  the  entire 

scatterer  and  denotes  the  average  energy  backscattered  from  each  of 

the  M  subscatterers.  Since  the  energy  distribution  backscattered  across 

the  scatterer  is  assumed  uniformly  distributed,  the  relationship 

E„  *  E/M  holds. 

M 

V.  DETECTION  PROBABILITY  -  TARGET  OF  M  RANGE  CELLS 

A.  SINGLE  PULSE  WAVEFORM  (N  -  1) 

First  we  consider  a  waveform  consisting  of  one  single  pulse.  Three 
target  models  are  evaluated  having  the  following  subscatterers. 

1)  Constant  amplitude  scatterers 

2)  Rayleigh  scatterers 

3)  Dominant  plus  Rayleigh  scatterers 

In  the  following  discussion,  the  probability  of  detection  is 
plotted  as  a  function  of  the  signal-to-noise  ratio  energy  per  pulse 
backscattered  from  the  target. 

1)  Constant  Amplitude  Scatterers 

A  scatterer  which  has  a  nonrandom  amplitude  and  a  uniformly 
distributed  phase  is  denoted  a  constant  amplitude  scatterer.  For 
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constant  amplitude  scatterers,  the  detection  probability,  P^,  for  at 

least  one  hit  out  of  M  hits,  each  hit  having  a  detection  probability 

of  is  given  by  equation  (3)  where 

DM 

PDM  =  1  "T\JV[  °>  (5)’ 

E/M  is  the  average  energy  in  the  received  signal  from  each  subscatterer 

and  Tylijjjk1]  is  the  incomplete  Toronto  function  [3],  [A],  [18]  -  [20] 

2 

with  threshold  to  noise  energy  ratio,  y  ,  and  arguments  i,  j,  and  k. 

For  the  constant  amplitude  case,  E  =  E.  A  four  term  expansion  of  the 
incomplete  Toronto  function  is  carried  out  in  appendix  A.  Appendix  B 
provides  a  discussion  on  how  the  incomplete  Toronto  function  relates  to 
probability  theory. 

For  constant  amplitude  subscatterers,  the  probability  of  detection 
as  a  function  of  the  signal-to-noise  ratio  is  shown  in  figures  3,  A, 
and  5.  For  each  value  of  M,  the  constant  amplitude  scatterer  model  has 
been  reconfigured  into  M  constant  amplitude  subscatterers  each  randomly 
located  within  each  of  the  M  range  cells.  Thus,  each  of  the  curves  in 
these  figures  represents  a  different  structure,  e.g.  one  sphere  for 
M  =  1,  two  spheres  for  M  =  2,  etc.  The  only  common  characteristics  are 
false  alarm  probability,  P^,  and  total  signal  energy,  E.  The  detection 
probability  degrades  as  the  resolution  index,  M,  increases  for  all 
values  of  signal-to-noise  ratio.  This  results  in  an  apparent 


signal-to-noise  ratio  loss 


2)  Rayleigh  Scatterers 


For  Rayleigh  subscatterers,  the  detection  probability,  P^,  for  at 
least  one  hit  out  of  M  possible  hits,  each  hit  having  a  detection 
probability  of  P^,  is  given  by  equation  (3)  where 


PDM  "  exp[  ‘ 


(  1  +  E/2MNq) 


For  Rayleigh  subscatterers  the  probability  of  detection  as  a 
function  of  the  signal-to-noise  ratio  is  shown  in  figures  6,  7,  and  8. 
Since  for  each  value  of  M  the  scatterer  model  is  complex  Gaussian 
(Rayleigh),  then  each  of  the  curves  in  these  figures  may  represent  the 
same  scatterer  by  virtue  of  the  property:  the  superposition  of  Gaussian 
random  variables  results  in  a  Gaussian  random  variable.  It  can  be  seen 
that  an  improvement  in  detection  probability  results  as  the  resolution 
index,  M,  increases  for  larger  values  of  signal-to-noise  ratio.  An 
increased  resolution  index,  M,  results  in  an  increase  in  gain  for  larger 
values  of  signal-to-noise  ratio  but  results  in  a  loss  for  the  lower 
signal-to-noise  ratios. 


3)  Dominant  Plus  Rayleigh  Scatterers 


For  the  dominant  plus  Rayleigh  subscatterers,  the  detection 
probability,  P^,  for  at  least  one  hit  out  of  M  possible  hits,  each  hit 
having  a  detection  probability  of  P^,  is  given  by  equation  (3)  where 
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For  dominant  plus  Rayleigh  subscatterers  the  probability  of 
detection  as  a  function  of  the  signal-to-noise  ratio  is  shown  in 
figures  9,  10,  and  11.  Again  the  scatterer  has  been  reconfigured  into  M 
dominant  plus  Rayleigh  subscatterers  each  randomly  located  within  each 
of  the  M  range  cells.  M  dominant  plus  Rayleigh  subscatterers  do  not 
quite  superimpose  to  result  in  a  M  ■  1  dominant  plus  Rayleigh  target. 
Instead  a  Rayleigh  target  would  result.  These  curves,  however,  do 
represent  the  same  overall  false  alarm  probability,  P^a,  and  total 

signal  energy,  E,  over  the  length  of  the  target.  The  appropriate 
comparisons  are  discussed  in  the  next  paragraph.  In  these  figures  there 
is  an  apparent  slight  improvement  in  detection  probability  as  the 
resolution  index,  M,  increases  for  larger  values  of  signal-to-noise 
ratio.  For  lower  false  alarm  probabilities  a  larger  signal-to-noise 
ratio,  E/Nq,  is  required  in  order  to  maintain  the  same  detection 
probability.  This  scatterer  represents  a  case  where  the  dominant  part 
of  the  scatterer  acts  as  the  constant  amplitude  scatterer  of  1)  above 
and  the  "plus  Rayleigh"  part  acts  as  the  Rayleigh  scatterer  of  2)  above. 


4)  Realistic  Cases 

In  reality,  a  scatterer  consisting  of  equal  amplitude  and/or 
dominant  plus  Rayleigh  subscatterers  approaches  a  complex  Gaussian 
(Rayleigh)  scatterer  as  the  number  of  subscatterers  increases.  This  is 
validated  by  virtue  of  the  central  limit  theorem  and  holds  for  all  the 
subscatterer  models  considered  above.  To  evaluate  these  realistic  cases 
we  plot  the  detection  probability  of  the  Rayleigh  scatterer  for  M  *  1 
with  the  detection  probability  of  its  constituent  constant  amplitude 
subscatterers  for  M  *  5  and  10.  The  result  is  shown  in  figures  12,  13, 
and  14.  It  can  be  seen  that  for  resolution  indices  of  around  M  =  5,  a 
marked  improvement  in  performance  occurs  as  the  signal-to-noise  ratios 
increase  beyond  that  corresponding  to  80%  detection  probability.  It 
should  be  kept  in  mind  that  correspondingly,  a  loss  occurs  below  the  80% 
probability,  i.e.  at  low  signal-to-noise  ratios.  Similar  improvements 
for  the  Rayleigh  scatterer  consisting  of  dominant  plus  Rayleigh 
subscatterers  are  illustrated  in  figures  15,  16,  and  17. 


B.  MULTIPLE  ENVELOPE  PULSE  INTEGRATION  -  N  PULSES 

We  now  consider  waveforms  consisting  of  a  train  of  N  pulses  and  the 
scatterer  divided  into  M  subscatterers.  The  same  subscatterers  treated 
previously  will  be  considered  with  the  addition  of  their  fluctuation 
characteristics . 
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1)  Constant  amplitude  scatterers 

2)  Rayleigh  scatterers 


a)  Slow  fluctuating  target 

b)  Fast  fluctuating  target 

3)  Dominant  plus  Rayleigh  scatterers 

a)  Slow  fluctuating  target 

b)  Fast  fluctuating  target 

In  the  probability  of  detection  plots  which  follow,  the  independent 
variable  is  signal-to-noise  ratio  energy  per  pulse  backscattered  from 
the  target. 

1)  Constant  Amplitude  Scatterers 

The  detection  probability  for  the  constant  subscatterers  is  given 
by  equation  (3)  where 

- 1 

NE/MNq  ]  (8). 

For  the  constant  amplitude  case,  E  ■  E. 

For  the  constant  amplitude  target,  the  probability  of  detection  as 
a  function  of  the  signal-to-noise  ratio  is  shown  in  figures  18,  19, 
and  20  for  a  four  pulse  train,  N  ■  4,  and  in  figures  21,  22,  and  23  for 
N  =  16.  Again,  as  in  figures  3,  4,  and  5,  for  each  value  of  M,  the 
constant  amplitude  scatterer  model  has  been  reconfigured  into  M  constant 
amplitude  subscatterers  each  randomly  located  within  each  of  the  M  range 
cells.  Thus,  each  of  the  curves  in  these  figures  represents  a  different 
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scatterer  where  the  only  common  characteristics  are  false  alarm 
probability,  P^a,  total  signal  energy,  E,  and  the  number  of  pulses,  N. 
The  best  performance  is  achieved  for  a  resolution  index  of  H  =>  1. 
Detection  probability  decreases  for  any  signal-to-noise  ratio  as  the 
resolution  index,  M,  increases. 


2)  Rayleigh  Scatterers 


For  M  Rayleigh  subscatterers  the  detection  probability  is  given  by 
equation  (3)  where 

a)  for  a  slow  fluctuating  target 

PDM  "  1  -  if  »  N"21  +  f"1  +  - -  1  ><  expf-  - — -  *"] 

LnJn-T1  J  L  NE/2MNo  J  L  1  +  NE/2MNo  J 

Xif-— - - - - 3-,  N.21  (9) 
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b)  and  for  a  fast  fluctuating  target 
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![.,.]  is  the  incomplete  gamma  function. 
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For  the  Rayleigh  slow  fluctuating  target  the  probability  of 
detection  as  a  function  of  the  signal-to-noise  ratio  is  shown  in 
figures  24,  25,  and  26  for  N  =  4  and  in  figures  27,  28,  and  29  for 
N  =  16.  As  discussed  after  equation  (6),  each  of  the  curves  in  each 
figure  represents  detection  performance  derived  from  the  same  scatterer. 
Such  a  relationship  is  a  consequence  of  the  fact  that  the  sum  of 
Gaussians  is  a  Gaussian.  Note  that  in  figure  24  the  detection 
probability  increases  at  a  larger  rate  with  respect  to  signal-to-noise 
ratio  for  increasingly  higher  values  of  M.  In  figure  24  for  M  ■  2, 
there  is  improved  detection  performance  with  respect  to  M  ■  1  when  the 
signal-to-noise  ratio  exceeds  8.9  dB.  For  signal-tc-naise  ratios  above 
13.1  dB,  maximum  performance  is  obtained  for  a  resolution  index  of 
M  *  5.  Above  17.1  dB  maximum  performance  is  obtained  for  a  resolution 
index  of  M  *  10.  Similar  effects  are  noted  in  figures  25  through  29. 

For  lower  false  alarm  probabilities  these  crossover  points  occur  at  a 
higher  signal-to-noise  ratio.  On  the  other  hand,  detection  probability 
deteriorates  for  increasingly  higher  values  of  M  when  the 
signal-to-noise  ratio  is  low. 

For  the  Rayleigh  slow  fluctuating  target,  the  N  pulse  samples  from 
each  subscatterer  are  dependent;  therefore,  the  scatterer  amplitude 
estimate  at  the  output  of  the  integrator  has  the  same  proportionate 
variance  as  that  of  the  single  pulse  case  for  all  M.  Its  detection 
performance  has  a  similar  relationship  with  respect  to  M  as  that  of  the 
Rayleigh  target  single  pulse  case.  Compare  figures  24  through  29  with 
figures  6  through  8. 
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For  the  Rayleigh  fast  fluctuating  target  the  probability  of 
detection  as  a  function  of  the  signal-to-noise  ratio  is  shown  in 
figures  30,  31,  and  32  for  N  =  4  and  in  figures  33,  34,  and  35  for 
N  =  16.  Here,  the  amplitude  of  each  of  the  M  subscatterers  has  been 
estimated  by  inputting  N  independent  samples  into  the  N  pulse  integrator 
shown  in  figure  2.  As  N  increases  the  estimate,  Y(j),  approaches  that 
of  the  M  constant  amplitude  subscatterers.  Under  these  conditions,  the 
curves  for  each  M  represent  the  same  target.  As  the  number  of  pulses, 

N,  increases,  the  detection  probability  tends  to  degrade  for  all  values 
of  signal-to-noise  ratio  for  increasingly  higher  values  of  M  as  was  the 
case  for  the  constant  amplitude  scatterers  of  figures  18  through  23. 

Note,  the  improvement  that  was  derived  from  increasing  M  in  the 
Rayleigh  single  pulse  case  or  the  Rayleigh  slow  fluctuating  target  case 
does  not  occur  with  increasing  N  in  the  fast  fluctuating  target  case. 

For  example,  in  the  fast  fluctuating  scatterer  of  figure  30,  the  M  =  2 
case  outperforms  the  M  =  1  case  when  the  signal-to-noise  ratio  exceeds 
12.7  dB.  From  figure  24,  i.e.  the  Rayleigh  slow  fluctuating  scatterer, 
the  crossover  point  was  at  8.9  dB.  As  the  number  of  independent 
samples,  N,  increases  as  in  figure  33  through  35,  the  estimated 
amplitude  approaches  that  of  the  constant  amplitude  case,  i.e.  the 
error  variance  decreases.  Also,  the  detection  performance  approaches 
that  of  the  constant  amplitude  case.  Comparing  figures  30  through  35 
with  figures  18  through  23  a  striking  resemblance  is  observed  between 
the  N  pulse  fast  fluctuating  target  and  the  constant  amplitude  target. 
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3)  Dominant  Plus  Rayleigh  Scatterers 


where 
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For  the  slow  fluctuating  dominant  plus  Rayleigh  target  the 
probability  of  detection  as  a  function  of  the  signal-to-noise  ratio  is 
shown  in  figures  36,  37,  and  38  for  N  -  4  and  in  figures  39,  40,  and  41 
for  N  =  16.  The  different  curves  for  each  M  do  not  represent  the  same 
scatterer  as  in  the  previously  discussed  Rayleigh  case.  However  they  do 
represent  the  same  overall  false  alarm  probability,  ,  total  signal 
energy,  E,  and  number  of  pulses  integrated,  N.  Since  such  a  target 


consists  of  a  constant  amplitude  and  a  Rayleigh  component,  the  detection 
performance  falls  somewhere  between  that  of  the  constant  amplitude 
scatterer  of  equation  (8)  and  the  Rayleigh  slow  fluctuating  scatterer  of 
equation  (9).  In  figure  36  it  is  noted  that  the  detection  probability 
increases  at  a  higher  rate  as  the  resolution  index,  M,  increases.  Above 
the  signal-to-noise  ratio  of  11.1  dB,  maximum  performance  is  obtained 
for  a  resolution  index  of  M  =  2  while  above  15.0  dB  maximum  performance 
is  obtained  for  M  =  5.  A  deterioration  in  performance  is  noted  with 
increasing  M  for  low  values  of  signal-to-noise  ratio. 

For  the  fast  fluctuating  dominant  plus  Rayleigh  target  the 
probability  of  detection  as  a  function  of  the  signal-to-noise  ratio  is 
shown  in  figures  42,  43,  and  44  for  N  «  4  and  in  figures  45,  46,  and  47 
for  N  =  16,  Again  the  N  independent  pulse  samples  result  in  M  constant 
amplitude  estimates  at  the  output  of  the  N  pulse  integrator.  Therefore, 
the  curves  for  each  M  represent  the  same  target.  In  figure  42,  where 
N  =  4,  the  detection  probability  incurs  only  a  slight  increase  as  the 
resolution  index,  M,  increases.  Above  the  signal-to-noise  ratio  of 
14.2  dB,  maximum  signal-to-noise  ratio  is  obtained  for  a  resolution 
index  of  M  =  2.  In  figure  45,  where  N  =  16,  there  is  a  performance  loss 
for  all  signal-to-noise  ratios  with  an  increase  of  M. 

Note,  however,  that  the  same  relationships  as  with  the  Rayleigh 
targets  hold  although  not  as  pronounced.  The  integrated  N  dependent 
pulses  from  the  slow  fluctuating  target  has  a  higher  rate  of  improvement 
with  signal-to-noise  ratio  as  M  increases.  Integration  of  the  N 
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independent  pulses  from  a  fast  fluctuating  target  results  in  an 
effective  performance  loss  as  M  increases.  Here,  as  in  the  Rayleigh 
fast  fluctuating  case,  the  variance  of  the  amplitude  estimate  decreases 
with  an  increase  in  N,  the  number  of  independent  samples.  Thus  for 
increasing  N,  the  variable  Y(j)  and  performance  thereof  approaches  that 
of  the  constant  amplitude  scatterer  when  the  target  is  fast  fluctuating. 

A)  Realistic  Cases 

Detection  probability  curves  are  shown  in  figures  48  to  50  for 
number  of  pulses,  N  »  4,  and  in  figures  51  to  53  for  number  of  pulses 
N  =  16.  Performance  for  the  fast  and  slow  Rayleigh  fluctuating  targets, 
M  »  1,  is  compared  with  the  performance  of  constant  amplitude 
subscatterers,  M  ■  5  and  M  »  10.  One  can  observe  the  improvement/loss 
in  performance  obtained  when  either  a  Rayleigh  fast  or  slow  fluctuating 
target  is  resolved  into  M  *  5  and  10  constant  scatterer  range  cells. 

When  this  is  indeed  the  scenario,  the  Rayleigh  slow  fluctuating  target 
incurs  a  signal-to-noise  improvement  as  the  resolution  index,  M, 
increases,  the  number  of  pulses,  N,  increases,  and  the  signal-to-noise 
ratio  is  such  so  that  the  probability  of  detection  is  above  80%.  On  the 
other  hand,  as  N  increases,  a  fast  fluctuating  target  is  best  detected 
at  a  resolution  index  of  M  *  1. 

Figures  5A  through  59  show  detection  probability  versus 
signal-to-noise  ratio  for  slow  fluctuating  Rayleigh  targets,  i.e. 
integration  of  dependent  samples.  When  the  slow  fluctuating  Rayleigh 


-  36  - 


Probability  of  Detection  vs.  Signal  to  Noise  Ratio 


/S  /JO 


_ Figure  63 

—  Rayleigh  Scatterers 
_  Fast  M  ■  1 

—  Dominant  plus  Rayleigh  Scatterers 
— Fast  M  •  2,  5.  10 

_  K  ■  16  pulses 


Probability  of  Detection  vs.  Signal  to  Noise  Ratio 


i  ft  /io 


Figure  64 

Rayleigh  ScatUrers 
Fast  H  ■  1 

Doeinant  plus  Raylet*  ScatUrers 
Fast  N  •  2.  S,  16 

u  ■  16  pulses 


Signal  to  Noise  Ratio,  E/N  in  dfl 


Probability  of  Detection  vs.  Signal  to  Noise  Ratio 


Figure  66 

Rayleigh  ScatUrers 
Fast  *  ■  1 

Dominant  plus  Rayleigh  ScatUrers 
Fast  H  ■  2,  5,  10 

N  •  16  pulses 

?u  ■  1  X  10'8 


\r  V*.  *  *  * 


42 


For  high  signal-to-noi'se  ratios: 

o  Fast  fluctuating  Rayleigh  targets  used  in  conjunction  with 
large  N  pulse  waveforms  incur  a  performance  degradation  as  the 
number  of  range  cells,  M,  increases  to  the  point  at  which  the 
individual  subscatterers  are  resolved. 

o  Single  pulse  fast  fluctuating  Rayleigh  targets,  in  conjunction 
with  small-N  pulse  waveforms,  and  slow  fluctuating  Rayleigh 
targets,  with  any  N,  incur  a  performance  improvement  as  the 
number  of  range  cells,  M,  increases  to  the  point  at  which  the 
individual  subscatterers  are  resolved. 

For  low  signal-to-noise  ratios: 

o  The  best  aetection  performance  is  achieved  at  M  =  1,  i.e.  the 
entire  target  occupies  one  single  range  resolution  cell.  This 
applies  for  any  scatterer  model. 

When  these  procedures  are  applied  to  distributed  clutter,  the  total 
false  alarm  count  may  increase  or  decrease  depending  on  the  statistical 
properties  of  the  clutter.  Just  as  the  Rayleigh  slow  fluctuating  target 
incurs  an  increase  in  detection  probability  with  increased  range 
resolution,  so  would  Rayleigh  slow  fluctuating  clutter  incur  a  higher 
false  alarm  count  (given  a  similar  threshold  adjustment  for  constant 
false  alarm  due  to  thermal  noise)  over  the  same  range  span. 


APPENDIX  A 


Development  of  the  Gram-Charlier  Series  for  Detection 
Probability  of  Constant  Amplitude  Scatterers 


■A  search  was  made  to  determine  analytical  expressions/expansions 
for  the  detection  probability  of  constant  targets  in  noise;  the 
incomplete  Toronto  function  or  as  sometimes  denoted  Marcum's  Q  function. 
No  suitable  expansion  was  found  that  could  be  readily  programmed.  It 
was  decided  to  undertake  a  fourth-order  Gram-Charlier  series 
approximation.  The  data  generated  from  these  approximations  was 
verified  and  found  to  be  well  within  the  accuracy  of  the  curves 
appearing  in  published  literature. 

The  Gram-Charlier  series  expansion  is  given  by 
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and  where  p(x)  is  the  approximated  probability  density  function.  Using 
the  procedure  given  in  [5]  we  may  determine  the  coefficients  c^  as 
follows. 
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First  we  determine  the  moments  from  the  characteristic  function 
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where  m  is  the  n  th  moment  of  the  distribution  of  Y  and  C„(0)  is 
n  Y 

the  n  th  derivative  of  the  characteristic  function  C^( £ )  evaluated 
at  £  ■  0.  For  the  particular  case  of  interest,  i.e.  the  constant 
scatterer  in  Gaussian  noise,  the  characteristic  function  jR^)  of  Y 
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Rp  is  equal  co  2E/MNq.  In  the  constant  amplitude  case,  E 
first  derivative  of  C^( ^ ;Rp)  is 


=  E.  The 


d  cY(  £  J*  ) 


jNR  j  £  NR 

- P '  ■ —2  exp(  - ^  )  + 


•  b  \N+1 


j^NR 

exp(  - %  ) 

2(l-j A ) 
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,  v .  *  v 

j  ;-w.f 


Evaluating  m^  yields 


m  =  Y  =  (N  R  /2)  +  N 
1  P 


In  a  similar  manner  the  second  and  third  moments  are: 


m2  =  (N  Rp  /2)  +  2(N  Rp  /2)(W-1)  +  N(N+1) 


m3  =  (N  Rp  /2)3  +  3(N  Rp  /2)2(Nf2)  +  3(N  Rp  / 2) ( N+l ) ( N+2)  + 


+  N(N+l)(N+2) 


Coefficients  c^  can  be  evaluated  with  the  relation 


1  "  J 


p(x)  Hj[C  —r5  )  dx 


where  the  first  four  Hermite  polynomials  are 


Ho(x)  =  1 


HjCx)  =  x 


H2(x)  =  x  -  1 


H^Cx)  =  x  -  3x 


Evaluating  the  above  integral  yields 


c  =  1 
o 


Cj  -  c2  =  ° 


CL3 

C3  =  ~  ^  ~ 37  ^ 
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v  y  v\-v-v. .  *• 

V/'.v 


where  ctj  *s  the  *  th  central  moment  of  p(x)  normalized  by  cr1: 


Since 


then 


Cl 


M+OJ 

-  ~t  J  (x  - 
*  -00 

■/ 


p(x)  dx 


x  p(x)  dx 


CL3  ■  (  m3  "  3m2m1  +  2m^  )  /  cr3 


where 


2  2 
cr  -  m2  -  mj 


The  order  of  magnitude  of  coefficient  c^  does  not  decrease 
uniformly  as  i  increases.  The  Edgeworth  series  is  yielded  by  regrouping 
the  terms  of  the  Gram-Charlier  series  with  the  same  magnitude.  The 
hierarchy  of  terms  in  descending  order  of  magnitude  is 


i  =  0 
i  =  3 
i  *  4 ,  6 
i  -  5,  7,  9 


A  fourth-order  approximation  of  p(x)  is  generated  by  expanding  the 
first  four  terms  of  the  Gram-Charlier  series,  i.e.  i  =  0,  1,2,  3. 

Upon  substitution  of  all  the  above  terms  into  equation  (A-l)  we  have  for 
the  fourth-order  probability  density  function 


p(x)  = 


-  c^(x  -  3x)]  exp(  - 
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'  fc  v*  *  *  .  v-  \  v*  ■  w>  i v  v 

Vv>>V'V^J.S.-v'-V-V.'4 
W-  >  - -\-v V- y  •  \ 


yyryA  t. 


r 

♦ 

* 


k 

t 


The  detection  probability  is  obtained  from  the  integral 


J 


CO 


PDM=J  P(X)  dx 

x-x, 
b 


,  /• 00  3  2 

PDM  »  —  J  I=0  -  c3(x  ‘  3x»  exp(  -  T  >  dx 
^  *  xb=(Yb-Y)/<r 


pdm  *  I ec£c( V  ^ >  +  c3(y=r  > ( 1  -  > exp(  ■ 


2 


where 


I 


CO  2 

erfc(xb/  NjT1 )  =  2  ]  )  exp(  -  )  dx 


x=x,  \f 


2n 


is  the  complimentary  error  function.  The  error  function  is  readily 
evaluated  by  a  series  expansion.  (3) 


cO 


erf(z)  =  1  -  erfc(z)  =  -  2 

\F  n-0 


,  ,Nn  2n+l 

(-1)  z 
n !  ( 2n+l ) 


Also  erf(-z)  =  -  erf(+z) 


(3) 

See  this  reference  for  information  about  the  error  function. 
"Handbook  of  Mathematical  Functions,"  U.S.  Department  of 
Commerce,  National  Bureau  of  Standards,  AMS  55,  June  1964. 


APPENDIX  B 


The  Incomplete  Toronto  Function 


The  integral 


rB 

TB(m,n,r)  -  2  r"  m+1  exp(  -r2  ) J  tm  n  exp(  -t2  )  In(  2rt  )  dt 

*  0 

(  B— 1 

is  called  the  "Incomplete  Toronto  Function".  In  the  above  integral 

I  (x)  denotes  the  modified  Bessel  function  of  the  first  kind  and  of 
nN 

order  n. 


An  application  of  this  function  is  in  the  expression  for  the 
probability  of  declaring  that  the  signal  s^  is  present  (H^  hypothesis) 
The  observation  under  the  two  hypothesis  is 


:  H 


si  +  ni 


:  H, 


for  i  =  1,  2,  3,  . N >  The  n^  are  complex  Gaussian,  zero  mean,  and 
variance  of  Nq.  The  s^  are  constant  amplitude,  E,  with  uniform  phase; 


-  IT  <  9.  <  +  TT  . 
a 


A  likelihood  ratio  test  for  the  above  hypothesis  leads  to  the 


following  decision  rule: 


where 


N 

Y  ■  5  2  IVjl2  (B-2) 

i-1 

and  is  a  threshold  which  establishes  the  boundary  and  the  critical 
region. 

The  probability  density  function  of  Y  given  E  and  given  Nq  is 

p(Y:Rp)  »  C2Y/NRp)(N_1)/2  exp(-Y-NRp/2)  IN_1(\|  2NRpY  ^  5  Y  >  0 

*  0  :  Otherwise. 

(B-3) 

where  is  the  modified  Bessel  function  of  the  first  kind,  of 

order  N-l,  and  Rp  =»  2E/Nq.  (4)  For  the  constant  amplitude  case, 

E  =  E.  The  probability  that  Y  _>  Y^,  i.e.  the  critical  region,  is  the 
cumulative  density  function 


*  CO 

Pd  =J  (2Y/NRp)(N"1)/2  exp(-Y  -  NRp/2)  2NRpY 1 )  dY 


(B-4) 


(4) 

For  simplicity  we  confine  the  discussion  to  a  single  scatterer  and 
avoid  any  association  with  over  resolved  targets,  thus  reference  to 
resolution  index,  M,  has  been  omitted. 


This  integral  is  of  the  form  of  equation  (B-l).  We  may  express  IV 
in  terms  of  the  incomplete  Toronto  function. 


Pd  -  1  -  T^p-nC  2N-1,  N-l,  \[nT7T) 


(B-S) 


This  is  the  probability  that  the  Hq  hypothesis  will  be  rejected  or 
equivalently,  that  the  constant  amplitude  signal  will  be  detected. 

An  approximation  to  equation  (B-5)  may  be  obtained  using  a 
Gram-Charlier  series.  A  four  term  approximation  is  as  follows. 


Pd  “  \  erfc<V  Nj7)  +  c3(  )  (  1  -  \2>  exP(  -  \  ) 


where 


erfc(x^/ 


\F>  ■ 2 J ” 

*  X“X, 


(  )  exp(  -  ) 


NF1 


dx 


is  the  complimentary  error  function. 
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